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INTRODUCTION TO B PHYSICS 



MATTHIAS NEUBERT 
Newman Laboratory of Nuclear Studies 
Cornell University, Ithaca, NY 14853, USA 



These lectures provide an introduction to various topics in heavy-flavor physics. 
We review the theory and phenomenology of heavy-quark symmetry, exclusive 
weak decays of B mesons, inclusive decay rates, and some rare B decays. 

1 Introduction 

The rich phenomenology of weak decays has always been a source of informa- 
tion about the nature of elementary particle interactions. A long time ago, /3- 
and /x-decay experiments revealed the structure of the effective flavor-changing 
interactions at low momentum transfer. Today, weak decays of hadrons con- 
taining heavy quarks arc employed for tests of the Standard Model and mea- 
surements of its parameters. In particular, they offer the most direct way 
to determine the weak mixing angles, to test the unitarity of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix, and to explore the physics of CP violation. 
Hopefully, this will provide some hints about New Physics beyond the Stan- 
dard Model. On the other hand, hadronic weak decays also serve as a probe of 
that part of strong-interaction phenomenology which is least understood: the 
confinement of quarks and gluons inside hadrons. 

The structure of weak interactions in the Standard Model is rather simple. 
Flavor-changing decays are mediated by the coupling of the charged current 
J^r^ to the VF-boson field: 



^cc 




(1) 



where 



JCC = {Ve-,V^,Vr)^^ 




+ (WL,CL,iL)7''^CKM Sl 




(2) 



contains the left-handed lepton and quark fields, and 
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is the CKM matrix. At low energies, the charged-current interaction gives rise 
to local four-fermion couplings of the form 



£eff - -2V2GF JnnJnn.u : 



where 



1.16639(2) GeV"^ 



(4) 
(5) 



is the Fermi constant. 
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Figure 1: Examples of leptonic (B — > t Vr), semi-leptonic (f?" — > D^e Ue), and non- 
Icptonic {B^^ — > D+tt") decays of B mesons. 

According to the structure of the charged-current interaction, weak decays 
of hadrons can be divided into three classes: leptonic decays, in which the 
quarks of the decaying hadron annihilate each other and only leptons appear 
in the final state; semi-leptonic decays, in which both leptons and hadrons 
appear in the final state; and non-leptonic decays, in which the final state 
consists of hadrons only. Representative examples of these three types of decays 
are shown in Fig. H. The simple quark-line graphs shown in this figure are a 
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gross oversimplification, however. In the real world, quarks are confined inside 
hadrons, bound by the exchange of soft gluons. The simplicity of the weak 
interactions is overshadowed by the complexity of the strong interactions. A 
complicated interplay between the weak and strong forces characterizes the 
phenomenology of hadronic weak decays. As an example, a more realistic 
picture of a non-leptonic decay is shown in Fig. 

d 




Figure 2: More realistic representation of a non-leptonic decay. 

The complexity of strong-interaction effects increases with the number of 
quarks appearing in the final state. Bound-state effects in leptonic decays 
can be lumped into a single parameter (a "decay constant"), while those in 
semi-leptonic decays are described by invariant form factors depending on the 
momentum transfer between the hadrons. Approximate symmetries of the 
strong interactions help us to constrain the properties of these form factors. 
Non-leptonic weak decays, on the other hand, are much more complicated to 
deal with theoretically. Only very recently reliable tools have been developed 
that allow us to control the complex QCD dynamics in many two-body B 
decays using a heavy-quark expansion. 

Over the last decade, a lot of information on heavy-quark decays has been 
collected in experiments at e^e~ storage rings operating at the T(4s) reso- 
nance, and more recently at high-energy e^e~ and hadron colliders. This has 
led to a rather detailed knowledge of the flavor sector of the Standard Model 
and many of the parameters associated with it. In the years ahead the B fac- 
tories at SLAG, KEK, Cornell, and DESY will continue to provide a wealth of 
new results, focusing primarily on studies of CP violation and rare decays. 
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The experimental progress in heavy-flavor physics has been accompanied 
by a significant progress in theory, which was related to the discovery of heavy- 
quark symmetry, the development of the heavy-quark effective theory, and 
more generally the establishment of various kinds of heavy-quark expansions. 
The excitement about these developments rests upon the fact that they allow 
model-independent predictions in an area in which "progress" in theory often 
meant nothing more than the construction of a new model, which could be 
used to estimate some strong-interaction hadronic matrix elements. In Sec. ^ 
we review the physical picture behind heavy-quark symmetry and discuss the 
construction, as well as simple applications, of the heavy-quark effective theory. 
Section || deals with applications of these concepts to exclusive weak decays of 
B mesons. Applications of the heavy-quark expansion to inclusive B decays 
are reviewed in Sec. ^. We then focus on the exciting field of rare hadronic B 
decays, concentrating on the example of the decays B ttK. In Sec. we 
discuss the theoretical description of these decays and explain various strategies 
for constraining and determining the weak, CP-violating phase 7 = arg(V^f,) of 
the CKM matrix. In Sec. we discuss how rare decays can be used to search 
for New Physics beyond the Standard Model. 



2 Heavy-Quark Symmetry 

This section provides an introduction ±jOL±he ideas of heavy-quark symmetry0^ 
and the heavy-quark effective theory lItIJ, which provide the modern theoret- 
ical framework for the description of the properties and decays of hadrons 
containing a heavy quark. For a more detailed description of this subject, the 
reader is referred to the review articles in Refs. 18-24. 



2.1 The Physical Picture 

There are several reasons why the strong interactions of hadrons containing 
heavy quarks are easier to understand than those of hadrons containing only 
light quarks. The first is asymptotic freedom, the fact that the effective cou- 
pling constant of QCD becomes weak in processes with a large_momentum 
transfer, corresponding to interactions at short distance scales Ej'O. At large 
distances, on the other hand, the coupling becomes strong, leading to non- 
perturbative phenomena such as the confinement of quarks and gluons on a 
length scale i?had 1/Aqcd ^ 1 fm, which determines the size of hadrons. 
Roughly speaking, Aqcd ^ 0.2 GeV is the energy scale that separates the 
regions of large and small coupling constant. When the mass of a quark Q 
is much larger than this scale, mg 3> Aqcd, it is called a heavy quark. The 
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quarks of the Standard Model fall naturally into two classes: up, down and 
strange are light quarks, whereas charm, bottom and top are heavy quarks :| 
For heavy quarks, the effective coupling constant asimq) is small, implying 
that on length scales comparable to the Compton wavelength Ag ~ 1/™q 
the strong interactions are perturbative and much like the electromagnetic 
interactions. In fact, the quarkonium systems (QQ), whose size is of order 
Xq /asimq) <^ -Rhad, are very much hydrogen-hke. 

Systems composed of a heavy quark and other light constituents are more 
complicated. The size of such systems is determined by i?had, and the typical 
momenta exchanged between the heavy and light constituents are of order 
■'^QCD- The heavy quark is surrounded by a complicated, strongly interacting 
cloud of light quarks, antiquarks and gluons. In this case it is the fact that 
Xq ^ -Rhad: i-e- that the Compton wavelength of the heavy quark is much 
smaller than the size of the hadron, which leads to simplifications. To resolve 
the quantum numbers of the heavy quark would require a hard probe; the soft 
gluons exchanged between the heavy quark and the light constituents can only 
resolve distances much larger than Xq. Therefore, the light degrees of freedom 
are blind to the flavor (mass) and spin orientation of the heavy quark. They 
experience only its color field, which extends over large distances because of 
confinement. In the rest frame of the heavy quark, it is in fact only the electric 
color field that is important; relativistic effects such as color magnetism vanish 
as rriQ — s- oo. Since the heavy-quark spin participates in interactions only 
through such relativistic effects, it decouples. 

It follows that, in the limit mg — > oo, hadronic systems which differ only 
in the flavor or spin quantum numbers of theJieavy quark have the same con- 
figuration of their light degrees of freedoniEl lI. Although this observation still 
does not allow us to calculate what this configuration is, it provides relations 
between the properties of such particles as the heavy mesons B, D, B* and D* , 
or the heavy baryons At, and Ac (to the extent that corrections to the infinite 
quark-mass limit are small in these systems). These relations result from some 
approximate symmetries of the effective strong interactions of heavy quarks 
at low energies. The configuration of light degrees of freedom in a hadron 
containing a single heavy quark with velocity v does not change if this quark 
is replaced by another heavy quark with different flavor or spin, but with the 
same velocity. Both heavy quarks lead to the same static color field. For Nh 
heavy-quark flavors, there is thus an ^\]{2Nh) spin-flavor symmetry group, un- 
der which the effective strong interactions are invariant. These symmetries are 
in close correspondence to familiar properties of atoms. The flavor symmetry 

"Ironically, the top quark is of no relevance to our discussion here, since it is too heavy 
to form hadronic bound states before it decays. 
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is analogous to the fact that different isotopes have the same chemistry, since 
to good approximation the wave function of the electrons is independent of 
the mass of the nucleus. The electrons only see the total nuclear charge. The 
spin symmetry is analogous to the fact that the hyperfine levels in atoms are 
nearly degenerate. The nuclear spin decouples in the limit me/rriN 0. 

Heavy-quark symmetry is an approximate symmetry, and corrections arise 
since the quark masses are not infinite. In many respects, it is complementary 
to chiral symmetry, which arises in the opposite limit of small quark masses. 
There is an important distinction, however. Whereas chiral symmetry is a 
symmetry of the QCD Lagrangian in the limit of vanishing quark masses, 
heavy-quark symmetry is not a symmetry of the Lagrangian (not even an ap- 
proximate one), but rather a symmetry of an effective theory that is a good 
approximation to QCD in a certain kinematic region. It is realized only in sys- 
tems in which a heavy quark interacts predominantly by the exchange of soft 
gluons. In such systems the heavy quark is almost on-shell; its momentum fluc- 
tuates around the mass shell by an amount of order Aqcd- The corresponding 
fluctuations in the velocity of the heavy quark vanish as Aqcd /"m-Q — *■ 0. The 
velocity pbecomes a conserved quantity and is no longer a dynamical degree of 
freedomllil. Nevertheless, results derived on the basis of heavy-quark symme- 
try are model-independent consequences of QCD in a well-defined limit. The 
symmetry-breaking corrections can be studied in a systematic way. To this 
end, it is however necessary to cast the QCD Lagrangian for a heavy quark, 

CQ^Q{tp-mQ)Q, (6) 

into a form suitable for taking the limit toq — s- oo. 

2.2 Heavy- Quark Effective Theory 

The effects of a very heavy particle often become irrelevant at low energies. It 
is then useful to construct a low-energy effective theory, in which this heavy 
particle no longer appears. Eventually, this effective theory will be easier to 
deal with than the full theory. A familiar example is Fermi's theory of the 
weak interactions. For the description of the weak decays of hadrons, the 
weak interactions can be approximated by point-like four-fermion couplings, 
governed by a dimensionful coupling constant Gp [cf. (jj)]. The effects of the 
intermediate W bosons can only be resolved at energies much larger than the 
hadron masses. 

The process of removing -the degrees of freedom of a heavy particle in- 
volves the following steps EZl CJ: one first identifies the heavy-particle fields 
and "integrates them out" in the generating functional of the Green functions 
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of the theory. This is possible since at low energies the heavy particle does not 
appear as an external state. However, whereas the action of the full theory 
is usually a local one, what results after this first step is a non-local effective 
action. The non-locality is related to the fact that in the full theory the heavy 
particle with mass M can appear in virtual processes and propagate over a 
short but finite distance Ax ~ 1/M. Thus, a second step is required to obtain 
a local effective Lagrangian: the non-local effective action is rewritten 
infinite series of local terms in an Operator Product Expansion (OPE) E3Ej. 
Roughly speaking, this corresponds to an expansion in powers of 1/M. ft is 
in this step that the short- and long-distance physics is disentangled. The 
long-distance physics corresponds to interactions at low energies and is the 
same in the full and the effective theory. But short-distance effects arising 
from quantum corrections involving large virtual momenta (of order M) are 
not described correctly in the effective theory once the heavy particle has been 
integrated out. In a third step, they have to be added in a perturbative way 
using renormalization-group techniques. These short-distance effects lead to 
a renormalization of the coefficients of the local operators in the effective La- 
grangian. An example is the effective Lagrangian for non-leptonic weak decays, 
in which radiative corrections from hard gluons with virtual momenta in the 
range between m\y and some low renormalization scale fi give rise, to. Wilson 
coefficients, which renormalize the local four-fcrmion interactions E^tJ. 

The heavy-quark effective theory (HQET) is constructed to provide a sim- 
plified description of processes where a heavy quark interacts^ith light degrees 
of freedom predominantly by the exchange of soft gluons 13 "liZl. Clearly, mg is 
the high-energy scale in this case, and Aqcd is the scale of the hadronic physics 
we are interested in. The situation is illustrated in Fig. ^. At short distances, 
i.e. for energy scales larger than the heavy-quark mass, the physics is pertur- 
bative and described by conventional QCD. For mass scales much below the 
heavy-quark mass, the physics is complicated and non-perturbative because 
of confinement. Our goal is to obtain a simplified description in this region 
using an effective field theory. To separate short- and long-distance effects, we 
introduce a separation scale /z such that Aqcd ^ M ^ ^^Q- The HQET will 
be constructed in such a way that it is equivalent to QCD in the long-distance 
region, i.e. for scales below /i. In the short-distance region, the effective the- 
ory is incomplete, since some high-momentum modes have been integrated out 
from the full theory. The fact that the physics must be independent of the 
arbitrary scale fj, allows us to derive renormalization-group equations, which 
can be employed to deal with the short-distance effects in an efficient way. 

Compared with most effective theories, in which the degrees of freedom 
of a heavy particle are removed completely from the low-energy theory, the 
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Figure 3: Philosophy of the heavy-quark effective theory. 

HQET is special in that its purpose is to describe the properties and decays of 
hadrons which do contain a heavy quark. Hence, it is not possible to remove 
the heavy quark completely from the effective theory. What is possible is to 
integrate out the "small components" in the full heavy-quark spinor, which 
describe the fluctuations around the mass shell. 

The starting point in the construction of the HQET is the observation that 
a heavy quark bound inside a hadron moves more or less with the hadron's 
velocity v and is almost on-shell. Its momentum can be written as 

= mQv'^ + fc" , (7) 

where the components of the so-called residual momentum k are much smaller 
than mq. Note that t; is a four- velocity, so that = 1. Interactions of the 
heavy quark with light degrees of freedom change the residual momentum by 
an amount of order Afc ~ Aqcd, but the corresponding changes in the heavy- 
quark velocity vanish as AqcD/mg — > 0. In this situation, it is appropriate to 
introduce large- and small-component fields, hy and Hy, by 

hy{x) = e'™«"-" P+ Q{x) , Hy{x) = e*"«''-" P_ Q{x) , (8) 



EW theory 
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where P+ and P_ are projection operators defined as 



P± 



1 ±^ 
2 



(9) 



It follows that 



Q{x) 



— imQV-x 



(10) 



Because of the projection operators, the new fields satisfy fh^ — hy and ^ — 
—Hy. In the rest frame, i.e. for = (1,0,0,0), hy corresponds to the upper 
two components of Q, while Hy corresponds to the lower ones. Whereas hy 
annihilates a heavy quark with velocity w, Hy creates a heavy antiquark with 
velocity v. 

In terms of the new fields, the QCD Lagrangian (^ for a heavy quark 
takes the form 

— hy iv-Dhy — Hy iiv-D + 2mQ) Hy + hy ip±Hy + Hy ip±hy , (11) 

where D'^ = Z?'' — v^^v-Dis orthogonal to the heavy-quark velocity: v ■ D± — 0. 
In the rest frame, Z?^ — (0, D ) contains the spatial components of the covariant 
derivative. From (^ij), it is apparent that hy describes massless degrees of 
freedom, whereas Hy corresponds to fluctuations with twice the heavy-quark 
mass. These are the heavy degrees of freedom that will be eliminated in the 
construction of the effective theory. The fields arc mixed by the presence of the 
third and fourth terms, which describe pair creation or annihilation of heavy 
quarks and antiquarks. As shown in the first diagram in Fig. ^ in a virtual 
process, a heavy quark propagating forward in time can turn into an antiquark 
propagating backward in time, and then turn back into a quark. The energy of 
the intermediate quantum state hhH is larger than the energy of the incoming 
heavy quark by at least 2mQ. Because of this large energy gap, the virtual 
quantum fluctuation can only propagate over a short distance Ax ~ I/uiq. 
On hadronic scales set by i?had = 1/Aqcd, the process essentially looks like a 
local interaction of the form 



where we have simply replaced the propagator for Hy by l/2mQ. A more 
correct treatment is to integrate out the small-component field Hy, thereby 
deriving a non-local effective action for the large-component field hy, which 
can then be expanded in terms of local operators. Before doing this, let us 
mention a second type of virtual corrections involving pair creation, namely 
heavy-quark loops. An example is shown in the second diagram in Fig. 0. 



hy ip± 



1 



ip±h.. 



(12) 



2mQ 
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Heavy-quark loops cannot be described in terms of the effective fields and 
Hy, since the quark velocities inside a loop are not conserved and are in no 
way related to hadron velocities. However, such short-distance processes are 
proportional to the small coupling constant as{mQ) and can be calculated in 
perturbation theory. They lead to corrections that are added onto the low- 
energy effective theory in the renormalization procedure. 




Figure 4: Virtual fluctuations involving pair creation of heavy quarks. Time flows to the 

right. 

On a classical level, the heavy degrees of freedom represented by Hy can 
be eliminated using the equation of motion. Taking the variation of the La- 
grangian with respect to the field Hy , we obtain 

{ivD + 2mQ)Hy^ipj_hy. (13) 

This equation can formally be solved to give 

HV = T^ ^--—Ip^hy, (14) 

showing that the small-component field Hy is indeed of order l/mq. We can 
now insert this solution into (^l]) to obtain the "non-local effective Lagrangian" 

£cS = hyivD hy + hyip± — — —ip±hy. (15) 

ZniQ + tv-D 

Clearly, the second term corresponds to the first class of virtual processes 
shown in Fig. ^. 

It is possible to derive this Lagrangian in a more elegant way by manipulat- 
ing th£.^generating functional for QCD Green functions containing heavy-quark 
fieldsll3. To this end, one starts from the field redefinition ( |lO| ) and couples 
the large-component fields hy to external sources py. Green functions with an 
arbitrary number of hy fields can be constructed by taking derivatives with 
respect to py. No sources are needed for the heavy degrees of freedom repre- 
sented by Hy . The functional integral over these fields is Gaussian and can be 
performed explicitly, leading to the effective action 



•S'cff = dxCcs~i\nA, (16) 
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with CcS as given in (ITs) . The appearance of the logarithm of the determinant 



A = exp Q Tr In [2mQ + iv-D - iri\^ 



(17) 



is a quantum effect not present in the classical derivation presented above. 
However, in this case the determinant can be regulated in a gauge-invariant 
way, and by choosiiLg_the gauge v ■ A = Q one can show that In A is just an 
irrelevant constant EZItj. 

Because of the phase factor in (|l^) , the x dependence of the effective heavy- 
quark field hv is weak. In momentum space, derivatives acting on hv produce 
powers of the residual momentum k, which is much smaller than mq. Hence, 
the non-local effective Lagrangian (^5|) allows for a derivative expansion: 

1 °° - / iv ■ D\ " 
CeS = Kiv-D + y^hyip±\-- ) ip±hy . (18) 



^0 ^ ^""0 - 

Taking into account that contains a P+ projection operator, and using the 
identity 



P+ip^ip^P+= 



P+ , (19) 



where i{D^^ , D^] = gs G^'^ is the gluon field-strength tensor, one finds that 



1 

2mQ ' AniQ 



C^ff = hy iv D hy + hy {iDi_f K + -j^ K (Jf,,. G"" hy + 0{l/ml) . (20) 



In the limit mg — > oo, only the first term remains: 

Coo ~ hy iv- D hy . (21) 

This is the effective Lagrangian of the HQET. It gives rise to the Feynman 
rules shown in Fig. ^. 

Let us take a moment to study the symmetries of this LagrangianEil Since 
there appear no Dirac matrices, interactions of the heavy quark with gluons 
leave its spin unchanged. Associated with this is an SU(2) symmetry group, 
under which Coo is invariant. The action of this symmetry on the heavy-quark 
fields becomes most transparent in the rest frame, where the generators S"* of 
SU(2) can be chosen as 

^'"Ko ' [S\S^]^^e^^'^S\ (22) 



11 



Figure 5: Feynman rules of the HQET and a are color indices). A heavy quark with 
velocity v is represented by a double line. The residual momentum k is defined in (m). 



Here cr* are the Pauli matrices. An infinitesimal SU(2) transformation /i„ —f 
(1 + ie- S ) hy leaves the Lagrangian invariant: 

SCoo — hy [iv-D, ie-S]hy — . (23) 

Another symmetry of the HQET arises since the mass of the heavy quark does 
not appear in the effective Lagrangian. For heavy quarks moving at the 
same velocity, eq. ( |2l] ) can be extended by writing 

Nh _ 

C^^J2~^liv-Dhl^. (24) 

1=1 

This is invariant under rotations in flavor space. When combined with the 
spin symmetry, the symmetry groiup-is promoted to SU(2A'^;i). This is the 
heavy-quark spin- flavor symmetry HEJ. Its physical content is that, in the 
limit mg oo, the strong interactions of a heavy quark become independent 
of its mass and spin. 

Consider now the operators appearing at order l/mg in the effective La- 
grangian (pO|). They are easiest to identify in the rest frame. The first operator, 

Okin = 7^ hy {iD^f hy -> hy {l D f hy , (25) 

2mQ 2mQ 

is the gauge-covariant extension of the kinetic energy arising from the residual 
motion of the heavy quark. The second operator is the non-Abelian analogue 
of the Pauli interaction, which describes the color-magnetic coupling of the 
heavy-quark spin to the giuon field: 

Cmae = hy tJui/ G'"' hy — > hy S-Bchy . (26) 

4mQ rriQ 
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Here S is the spin operator defined in (|2^), and Bl = —^t^^^G^^ are the 
components of the color- magnetic field. The chromo-magnetic interaction is a 
relativistic effect, which scales like l/mg. This is the origin of the heavy-quark 
spin symmetry. 



2.3 The Residual Mass Term and the Definition of the Heavy- Quark Mass 

The choice of the expansion parameter in the HQET, i.e. the definition of the 
heavy-quark mass toq, deserves some comments. In the derivation presented 
earlier in this section, we chose mq to be the "mass in the Lagrangian" , and 
using this parameter in the phase redefinition in (|l^ we obtained the effective 
Lagrangian (|2l|), in which the heavy-quark mass no longer appears. However, 
this treatment has its subtleties. The symmetries of the HQET allow a "resid- 
ual mass" 5m for the heavy quark, provided that 5m is of order Aqcd and 
is the same for all heavy-quark flavors. Even if we arrange that such a mass 
term is not present at the tree level, it will in general be induced by quantum 
corrections. (This is unavoidable if the theory is regulated with a dimensionful 
cutoff.) Therefore, instead of (|2l|) we should write the effective Lagrangian in 



the more general form 
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Coo = hyiv-Dhy — 5mhjjhy . (27) 

If we redefine the expansion parameter according to mg mq + Am, the 
residual mass changes in the opposite way: 5m 5m — Am. This implies that 
there is a unique choice of the expansion parameter mq such that 5m — 0. 
Requiring 5m = 0, as it is usually done implicitly in the HQET, defines-a-heavy- 
quark mass, which in perturbation theory coincides with the pole massE3. This, 
in turn, defines for each heavy hadron Hq a parameter A (sometimes called 
the "binding energy" ) through 

A = {mHQ - mq) . (28) 

mq — *oo 

If one prefers to work with another choice of the expansion parameter, the 
values of non-perturbative parameters such as A change, but at the same time 
one has to include the residual mass term in the HQET Lagrangian. It can 
be shown that the various parameters depending on the definition of mq en- 
ter the predictions for physical quantities in^iuch a way that the results are 
independent of the particular choice adopted t3. 

There is one more subtlety hidden in the above discussion. The quantities 
mq, A and 5m are non-perturbative parameters of the HQET, which have a 
similar status as the vacuum condensates in QCD phenomenology 123. These 
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parameters cannot be defined unambiguously in perturbation theory. The 
reason hes in the divergent behavior of perturbative expansions in large orders, 
which is associated with the existence oLaingularities along the real axis in the 
Borel plane, the so-called renormalons EJ^cZl. For instance, the perturbation 
series which relates the pole mass toq of a heavy quark to its bare mass, 

"^Q = "^Q '° (l + ci asimq) + C2 aKniQ) + . . . + c„ <(mQ) + . . . | , (29) 

contains numerical coefficients c„ that gnMz_.as n\ for large n, rendering the 
series divergent and not Borel summable The best one can achieve is to 
truncate the perturbation series at its minimal term, but this leads to an un- 
avoidable arbitrariness of order Atoq ~ Aqcd (the size of the minimal term) 
in the value of the pole mass. This observation, which at first sight seems a 
serious problem for QCD phenomenology, should not come as a surprise. We 
know that because of confinement quarks do not appear as physical states in 
nature. Hence, there is no unique way to define their on-shell properties such as 
a pole mass. Remarkably, QCD perturbation theory "knows" about its incom- 
pleteness and indicates, through the appearance of renormalon singularities, 
the presence of non-perturbative effects. One must first specify a scheme how 
to truncate the QCD perturbation series before non-perturbative statements 
such as Sm ~ become meaningful, and hence before non-perturbative pa- 
rameters such as mg and A become well-defined quantities. The actual values 
of these parameters will depend on this scheme. 

We stress that the "renormalon ambiguities" are not a conceptual problem 
for the heavy-quark expansion. In fact, it can be shown quite generally that 
these ambiguities cancel in all predictions for physical observablesEJ l3. The 
way the cancellations occur is intricate, however. The generic structure of the 
heavy-quark expansion for an observable is of the form: 

Observable - ClaJmn)] ( I + — + ...] , (30) 

\ J 

where C'lasirng)] represents a perturbative coefficient function, and A is a 
dimensionful non-perturbative parameter. The truncation of the perturba- 
tion series defining the coefficient function leads to an arbitrariness of order 
Aqcd/wq, which cancels against a corresponding arbitrariness of order Aqcd 
in the definition of the non-perturbative parameter A. 

The renormalon problem poses itself when one imagines to apply pertur- 
bation theory to very high orders. In practice, the perturbative coefficients are 
known to finite order in (typically to one- or two- loop accuracy) , and to be 
consistent one should use them in connection with the pole mass (and A etc.) 
defined to the same order. 
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2-4 Spectroscopic Implications 

The spin-flavor symmetry leads to many interesting relations between the prop- 
erties of hadrons containing a heavy ouatk. The most direct consequences 
concern the spectroscopy of such states BeJ. In the limit uiq — > oo, the spin 
of the heavy quark and the total angular momentum j of the light degrees 
of freedom are separately conserved by the strong interactions. Because of 
heavy-quark symmetry, the dynamics is independent of the spin and mass of 
the heavy quark. Hadronic states can thus be classified by the quantum num- 
bers (flavor, spin, parity, etc.) of their light degrees of freedomEJ. The spin 
symmetry predicts that, for fixed j ^ 0, there is a doublet of degenerate states 
with total spin J = j ± i . The flavor symmetry relates the properties of states 
with different heavy-quark flavor. 

In general, the mass of a hadron Hq containing a heavy quark Q obeys an 
expansion of the form 

"^//Q = "^Q + A + ^ + 0(l/m|) . (31) 

The parameter A represents contributions arising-from terms in the Lagrangian 
that are independent of the heavy-quark mass £3, whereas the quantity Am^ 
originates from the terms of order l/rUQ in the effective Lagrangian of the 
HQET. For the ground-state pseudoscalar and vector mesons, one can parame- 
trize the contributions from the kinetic energy and the chromo-rtLagnetic in- 
teraction in terms of two quantities Ai and A2, in such a way thatEj 



Am^ ^ -Ai + 2 



J(J+l)-§ A2. (32) 



The hadronic parameters A, Ai and A2 are independent of toq. They charac- 
terize the properties of the light constituents. 

Consider, as a flrst example, the SU(3) mass splittings for heavy mesons. 
The heavy-quark expansion predicts that 

ruBs - = As ~ + 0{l/mb) , 

rriDs - rriDa = A^ - A^ + 0{l/mc) , (33) 

where we have indicated that the value of the parameter A depends on the 
flavor of the light quark. Thus, to the extent that the charm and bottom 
quarks can both be considered sufficiently heavy, the mass splittings should be 
similar in the two systems. This prediction is confirmed experimentally, since 

mss -ruB^ = (90 ± 3) MeV , 

mDs - = (99 ± 1) MeV. (34) 
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As a second example, consider the spin splittings between the ground-state 
pseudoscalar (J — 0) and vector (J = 1) mesons, which are the members of 
the spin- doublet with j = |. From (|||) and (H), it follows that 

m%, - m% = 4A2 -I- 0{l/mb) , 

to|,. -m|j = 4A2 + 0(l/mc). (35) 
The data are compatible with this: 

m%. - m| w 0.49 GeV^ , 

ml, ~m\, K 0.55 GeV^ . (36) 

Assuming that the B system is close to the heavy-quark limit, we obtain the 
value 

A2 ~ 0.12 GeV^ (37) 

for one of the hadronic parameters in (p^). This quantity plays an important 
role in the phenomenology of inclusive decays of heavy hadrons. 

A third example is provided by the mass splittings between the ground- 
state mesons and baryons containing a heavy quark. The HQET predicts that 

"^Ai, - rnB = Abaryon - Amcson + 0{l/mb) , 

"^Ac - "Id = Abaryon - A,„cson + 0(1/771^) . (38) 

This is again consistent with the experimental results 
ruA, -ruB = (345 ± 9) MeV , 

TOA, - = (416± 1) MeV, (39) 

although in this case the data indicate sizeable symmetry-breaking corrections. 
The dominant correction to the relations ( p8| ) comes from the contribution of 
the chromo-magnetic interaction to the masses of the heavy mesons|| which 
adds a term 3A2/2toq on the right-hand side. Including this term, we obtain 
the refined prediction that the two quantities 

rriA, -me- ^ = (311 ± 9) MeV , 

niA -mo ^ = (320±l)MeV (40) 

2mD 



''Because of spin symmetry, there is no such contribution to the masses of Aq baryons. 
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should be close to each other. This is clearly satisfied by the data. 

The mass formula ( |3l| ) can also be used to derive information on the 
heavy-quark masses from the observed hadron masses. Introducing the "spin- 
averaged" meson masses ms — \{^b + SmB*) ~ 5.31 GeV and riiD = 
jimD + Smo' ) w 1.97 GeV, we find that 

TUb - TUc = {niB - mn) 1 1 - + 0{l/m%) \ . (41) 

1^ 2m,BmD J 

Using theoretical estimates for the parameter Ai, which lie in the range 

Ai = -(0.3±0.2) GeV^ (42) 

this relation leads to 



nib -mc== (3.39 ±0.03 ±0.03) GeV, (43) 

where the first error reflects the uncertainty in the value of Ai, and the second 
one takes into account unknown higher-order corrections. The fact that the 
difference of the pole masses, mb — rric, is known rather precisely is important 
for the analysis of inclusive decays of heavy hadrons. 



3 Exclusive Semi-Leptonic Decays 

Semi-leptonic decays of B mesons have received a lot of attention in recent 
years. The decay channel B D*£D has the largest branching fraction of 
all _B-meson decay modes. From a theoretical point of view, semi-leptonic 
decays are simple enough to allow for a reliable, quantitative description. The 
analysis of these decays provides much information about the strong forces 
that bind the quarks and gluons into hadrons. Schematically, a semi-leptonic 
decay process is shown in Fig. ^ The strength of the b c transition vertex 
is governed by the element Vcb of the CKM matrix. The parameters of this 
matrix are fundamental parameters of the Standard Model. A primary goal of 
the study of semi-leptonic decays of B mesons is to extract with high precision 
the values of \Vcb\ and \Vub\- We will now discuss the theoretical basis of such 
analyses. 

3.1 Weak Decay Form Factors 

Heavy-quark symmetry implies relations between the weak decay form factors 
of heavy mesons, which are of particular interest. These relations have been 
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D, D*, D* 



Figure 6: Semi-leptonic decays of B mesons. 

derivecLby Isgur and Wise i, generalizuig ideas developed by Nussinov and 
WetzelH, and by Voloshin and ShifmanQia. 

Consider the elastic scattering of a -B meson, B{v) B{v'), induced by a 
vector current coupled to the b quark. Before the action of the current, the light 
degrees of freedom inside the B meson orbit around the heavy quark, which 
acts as a static source of color. On average, the b quark and the B meson have 
the same velocity v. The action of the current is to replace instantaneously 
(at time t = to) the color source by one moving at a velocity u', as indicated 
in Fig. 1^. If V — v' , nothing happens; the light degrees of freedom do not 
realize that there was a current acting on the heavy quark. If the velocities are 
different, however, the light constituents suddenly find themselves interacting 
with a moving color source. Soft gluons have to be exchanged to rearrange 
them so as to form a B meson moving at velocity v' . This rearrangement leads 
to a form-factor suppression, reflecting the fact that, as the velocities become 
more and more different, the probability for an elastic transition decreases. 
The important observation is that, in the limit mf, oo, the form factor can 
only depend on the Lorentz boost ^ — v ■ v' connecting the rest frames of the 
initial- and final-state mesons. Thus, in this limit a dimensionless probability 
function ^(u • v') describes the transition. It is called the Isgur- Wise functionQ. 
In the HQET, which provides the appropriate framework for taking the limit 
rub oo, the hadronic matrix element describing the scattering process can 
thus be written as 

— {Biv')\K,j''b, \B{v)) = ■ v') [v + vY . (44) 
niB 

Here b^ and b^i are the velocity-dependent heavy-quark fields of the HQET. 
It is important that the function ^(v ■ v') does not depend on mf,. The factor 
1 / niB on the left-hand side compensates for a trivial dependence on the heavy- 
meson mass caused by the relativistic normalization of meson states, which is 
conventionally taken to be 

{Bip')\Bip))^2mBv'i2nf6^ip-p'). (45) 



18 



Note that there is no term proportional to [v — v')^ in (p^. This can be seen 
by contracting the matrix element with {v — w')^, which must give zero since 
fb^ — by and byif = byi. 




t<tQ t = tQ t>tQ 



Figure 7: Elastic transition induced by an external heavy-quark current. 

It is more conventional to write the above matrix element in terms of an 
elastic form factor Fc\{q'^) depending on the momentum transfer = (p — p')^- 

{B{v')\brb \B{v)) = Fei(q2) [p + p'r , (46) 

where p*^') = itibv^'-' ■ Comparing this with (^4|), we find that 

F,i{q^) = ■ v') , q^ = ~2mliv v' - I) . (47) 

Because of current conservation, the elastic form factor is normalized to unity 
at q^ = 0. This condition implies the normalization of the Isgur-Wise function 
at the kinematic point v ■ v' — 1, i.e. for v = v': 

e(l) = l. (48) 

It is in accordance with the intuitive argument that the probability for an 
elastic transition is unity if there is no velocity change. Since for v — v' the 
final-state meson is at rest in the rest frame of the initial meson, the point 
u • i>' = 1 is referred to as the zero-recoil limit. 

The heavy-quark flavor symmetry can be used to replace the b quark in 
the final-state meson by a c quark, thereby turning the B meson into a D 
meson. Then the scattering process turns into a weak decay process. In the 
infinite-mass limit, the replacement by/ Cy> is a symmetry transformation, 
under which the effective Lagrangian is invariant. Hence, the matrix element 

, ^ {D{v')\ c^n'^K \B{v)) - ^{v ■ v') {v + v'r (49) 

is still determined by the same function ^{v ■ v'). This is interesting, since 
in general the matrix element of a flavor-changing current between two pseu- 
doscalar mesons is described by two form factors; 

{D{v')\ crb \B{v)) - f+{q') (p + pT - /-(<?') {P - P'Y ■ (50) 
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Comparing the above two equations, we find that 



aw'), 



m\ + m?) — 2mBm£) v ■ v' . 



(51) 



Thus, the heavy-quark flavor symmetry relates two a priori independent form 
factors to one and the same function. Moreover, the normahzation of the 
Isgur-Wise function at v ■ v' = 1 now imphes a non-trivial normalization of 
the form factors f±{q^) at the point of maximum momentum transfer, (7^,ax = 



/±(9max) 



TUB ifTlB) 



(52) 



2y/mBmD 

The heavy-quark spin symmetry leads to additional relations among weak 
decay form factors. It can be used to relate matrix elements involving vector 
mesons to those involving pseudoscalar mesons. A vector meson with longi- 
tudinal polarization is related to a pseudoscalar meson by a rotation of the 
heavy-quark spin. Hence, the spin-symmetry transformation c^, c^,_relates 
B ^ D with B — > D* transitions. The result of this transformation isEl 



1 



y/mBrno' 
1 



{D*{v\e)\-c,,rK\B{v)) - ze^-^P elv'^vp ■ v') , 



^riiBniD* 



{D*{v',e)\c,,Yl5K\B{v)) = e*^ (w • w' + 1) V ^(ww'), 



(53) 



where e denotes the polarization vector of the D* meson. Once again, the 
matrix elements are completely described in terms of the Isgur-Wise function. 
Now this is even more remarkable, since in general four form factors, V{q^) for 
the vector current, and Ai{q^), i = 0, 1, 2, for the axial current, are required to 
parametetize these matrix elements. In the heavy-quark limit, they obey the 
relations eZI 



niB + rriD- 
2y/rnBrnD* 



1 - 



{rriB + moY 



Ai{q^). 



2mBm£it V ■ v' . 



(54) 
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Equations ( pT| ) and ( p4[ ) summarize the relations imposed by heavy-quark 
symmetry on the weak decay form factors describing the semi-leptonic decay 
processes B ^ Div and B D*iD. These relations are model-independent 
consequences of QCD in the limit where rrih, rric ^ Aqcd- They play a crucial 
role in the determination of the CKM matrix element \Vcb\- In terms of the 
recoil variable w = V^' , the differential semi-leptonic decay rates in the heavy- 
quark limit become E3 



Aw 



487r3 
G2 



dw 



■T^ \Vct\'' {rriB - niD^f ml, ^/w^ - I {w + if 



Aw rn\ ~ 2w ruBmu' -\- rn\,. 



eiw). (55) 



These expressions receive symmetry-breaking corrections, since the masses of 
the heavy quarks are not infinitely large. Perturbative corrections of order 
ck"(mQ) can be calculated order by order in perturbation theory. A more 
difhcult task is to control the non-perturbative power corrections of order 
(Aqcd/toq)"- The HQET provides a systematic framework for analyzing 
these corrections. For the case of weak-decay, form factors the analysis of 
the l/mg corrections was performed by Luketa. Later, Falk and the present 
author have analyzed the structure of 1/?71q corrections for both meson and 
baryon weak decay form factorsEl. We shall not discuss these rather technical 
issues in detail, but only mention the most important result of Luke's analy- 
sis. It concerns the zero-recoil limit, where an analogue of the Ademollo-Gatto 
theorem E2l can be proved. This is Luke's theorem!^, which states that the 
matrix elements describing the leading 1 /toq corrections to weak decay ampli- 
tudes vanisb at zero recoil. This theorem is valid to all orders in perturbation 
theory Most importantly, it protects the B-^ D*£D decay rate from 
receiving first-order 1 / mq corrections at zero recoilE3. [A similar statement is 
not true for the decay B ^ D £ D. The reason is simple but somewhat subtle. 
Luke's theorem protects only those form factors not multiplied by kinematic 
factors that vanish for v — v' . By angular momentum conservation, the two 
pseudoscalar mesons in the decay B D £9 must be in a relative p wave, and 
hence the amplitude is proportional to the velocity \vd\ of the D meson in the 
_B-meson rest frame. This leads to a factor (w^ — 1) in the decay ca±e. In such 
a situation, kinematically suppressed form factors can contribute EZI.] 
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3.2 Short- Distance Corrections 

In Sec. H, we have discussed the first two steps in the construction of the 
HQET. Integrating out the small components in the heavy-quark fields, a non- 
local effective action was derived, which was then expanded in a series of local 
operators. The effective Lagrangian obtained that way correctly reproduces 
the long-distance physics of the full theory (see Fig. ||). It does not contain 
the short-distance physics correctly, however. The reason is obvious: a heavy 
quark participates in strong interactions through its coupling to gluons. These 
gluons can be soft or hard, i.e. their virtual momenta can be small, of the order 
of the confinement scale, or large, of the order of the heavy-quark mass. But 
hard gluons can resolve the spin and flavor quantum numbers of a heavy quark. 
Their effects lead to a renormalization of the coefficients of the operators in 
the HQET. A new feature of such short-distance corrections is that through 
the running coupling constant they induce a logarithmic dependence on the 
heavy-quark massu. Since asimq) is small, these effects can be calculated in 
perturbation theory. 

Consider, as an example, the matrix elements of the vector current V — 
q^^Q. In QCD this current is partially conserved and needs no renormaliza- 
tion. Its matrix elements are free of ultraviolet divergences. Still, these matrix 
elements have a logarithmic dependence on toq from the exchange of hard 
gluons with virtual momenta of the order of the heavy-quark mass. If one goes 
over to the effective theory by taking the limit mq —f oo, these logarithms 
diverge. Consequeatly, the vector current in the effective theory does require 
a renormalizationtil. Its matrix elements depend on an arbitrary renormaliza- 
tion scale n, which separates the regions of short- and long-distance physics. 
If is chosen such that Aqcd ^ /i ^ mg, the effective coupling constant in 
the region between /x and mq is small, and perturbation theory can be used 
to compute the short-distance corrections. These corrections have to be added 
to the matrix elements of the effective theory, which contain the long-distance 
physics below the scale fi. Schematically, then, the relation between matrix 
elements in the full and in the effective theory is 

C (m ^) 

{V{mq))QCB = Coimq,fi) {Voiti))iiQET + — — — (Vi(^))hqet + . . . , 

(56) 

where we have indicated that matrix elements in the full theory depend on mq , 
whereas matrix elements in the effective theory are mass-independent, but do 
depend on the renormalization scale. The Wilson coefhcients Ci{mq, ^) are 
defined by this relation. Order by order in perturbation theory, they can 
be computed from a comparison of the matrix elements in the two theories. 
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Since the effective theory is constructed to reproduce correctly the low-energy 
behavior of the full theory, this "matching" procedure is independent of any 
long-distance physics, such as infrared singularities, non-perturbative effects, 
and the nature of the external states used in the matrix elements. 

The calculation of the coefficient functions in perturbation theory uses the 
powerful methods of the renormalization group. It is in principle straightfor- 
ward, yet in practice rather tedious. A comprehensive discussion of most of the 
existing calculations of short-distance corrections in the HQET can be found 
in Ref. 18. 



3.3 Model-Independent Determination of \Vc 



cb\ 



We will now discuss the most important application of the formalism de- 
scribed above in the context of semi-leptonic decays of B mesons. A model- 
independent determination of the CKM matrix element \Vcb\ based on heavy- 
quark symmetry can be obtained by laeasuring the recoil spectrum of D* 
mesons produced in i? — > D*iP decaysEa. In the heavy-quark limit, the dif- 
ferential decay rate for this process has been given in (|5^). In order to allow 
for corrections to that limit, we write 



dr 

dw 



487r3 



If 



1 



iKbP^'M, (57) 



where the hadronic form factor J-(w) coincides with the Isgur-Wise function 
up to symmetry-breaking corrections of order asimg) and A^qcd /itiq. The 
idea is to measure the product \Vcb \ ^{w) as a function of w, and to extract 
\Vcb\ from an extrapolation of the data to the zero- recoil point w — 1, where 
the B and the D* mesons have a common rest frame. At this kinematic point, 
heavy-quark symmetry helps us to calculate the normalization J-{1) with small 
and controlled theoretical errors. Since the range of w values accessible in this 
decay is rather small (1 < w < 1.5), the extrapolation can be done using an 
expansion around w = I: 



(58) 



The slope and the curvature c, and indeed more generally the complete 
shape of the .fn.r.nj factor, are tightly constrained by analyticity and unitarity 
requirements 0'L3. In the long run, the statistics of the experimental results 
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close to zero recoil will be such that these theoretical constraints will not be 
crucial to get a precision measurement of \Vcb\- They will, however, enable 
strong consistency checks. 



0.05 




Figure 8: CLEO data for the prtidiict \Vcb \ ^{w), as extracted from the recoil spectrum in 
B — ► D*£u decaysO- The Une shows a linear fit to the data. 



Measurements of the recoil spectrum have been performed by several ex- 
perimental groups. Figure^ shows, as an example, the data reported some time 
ago by the. CLEO Collaboration. The weighted average of the experimental 
results isEj 

= (35.2 ±2.6) X 10"^ (59) 

Heavy-quark symmetry implies that the general striicture of the symmetry- 
breaking corrections to the form factor at zero recoil is La 



^(1) = ??A 1 + Ox 



A 



QCD 



rriQ 



const X 



^^QCD 



. ^ = ?7A (1 + <5l/m2) : 



(60) 



where 77^ is a short-distance correction arising from the finite renormalization 
of the flavor-changing axial current at zero recoil, and ^i/m^ parameterizes 
second-order (and higher) power corrections. The absence of-first-order power 
corrections at zero recoil is a consequence of Lukels^theoremEJ. The one-loop 
expression for r]A has been known for a long timeHuQ: 



as{M) f nib + rric ^^rub 



\mb 



0.96. 



(61) 
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The scale M in the running couphng constant can be ftji^d by adopting the 
prescription of Brodsky, Lepage and Mackenzie (BLM) 13, where it is iden- 
tified with the average virtuality of the gluon in the one-loop diagrams that 
contribute to If as{M) is defined in the MS scheme, the result isEj M « 
0.51y/rri^mb. Several estimates of higher-order corrections to r]A have been 
discussed. A renormalization-group resummation of logarithms of tl 
{as In mh/mc)", as (as In m;,/mc)" and mc/mf,(as In mh/mc)" leads tot 
r]A ~ 0.985. On the other hand, a resummation of "renormalon-chain" con- 
tributions of the form where (3o = 11 — ^Uf is the first coefficient of 
the QCD /3-function, gives O rjA ~ 0.945. Using these partial resummations 
to estimate the uncertainty gives r/A — 0.965 ± 0.020. Recently, Czarnecki 
has improved this estimate by calculating rjA at two-loop ordered. His result, 
rjA — 0.960 ± 0.007, is in excellent agreement with the BLM-improved one- 



loop expression (61). Here the error is taken to be the size of the two-loop 
correction. 

The analysis of the power corrections is more difhcult, since it cannot 
rely on perturbation theory. Three approaches have been discussed: in the 
"exclusive approach" , all 1 /toq operators in the HQET are classified and their 
matrix elements estimated, leading to@@ <5i/m2 — — (3± 2)%; the "inclusive 
apprnach" has been used to derive the bound ^i/^a < —3%, and to estimate 
thatEZl (Si/to2 = —(7 ±3)%; the "hybrid approach" combines the virtues of thfi 
former two to obtain a more restrictive lower bound on Si/m^. This leads toH 
Si/m^ = -0.055 ±0.025. 

Combining the above results, adding the theoretical errors linearly to be 
conservative, gives 

J^(l) = 0.91 ±0.03 (62) 

for the normalization of the hadronic form factor at zero recoil. Thus, the 
corrections to the heavy-quark limit amount to a moderate decrease of the 
form factor of about 10%. This can be used to extract from the experimental 
result (^9|) the model-independent value 

\Vcb\ = (38.7 ± 2.8exp ± 1.3th) X 10-^ . (63) 



3.4 Measurements of B D*£ v and B —t D Form Factors and Tests of 
Heavy- Quark Symmetry 

We have discussed earlier in this section that heavy-quark symmetry implies 
relations between the semi-leptonic form factors of heavy mesons. They receive 
symmetry-breaking corrections, which can be estimated using the HQET. The 
extent to which these relations hold can be tested experimentally by comparing 
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the different form factors describing the decays B — ^ D^*'>£D at the same value 
of w. 

When the lepton mass is neglected, the differential decay distributions 
in B ^ D*£ V decays can be parameterized by three helicity amplitudes, or 
equivalently by three independent combinations of form factors. It has been 



suggested that a goocLch. 
the heavy-quark limit ESl 



lice for three such quantities should be inspired by 
1 One thus defines a form factor hAi{w), which up 
to symmetry-breaking corrections coincides with the Isgur-Wise function, and 
two form-factor ratios 



Ri{w) = 



R2{w) = 



1 - 



1 - 



{■niB + mD*Y 
{mB + rriD'Y 



Ai{q^)- 



(64) 



The relation between w and has been given in ( p^ ) . This definition is such 
that in the heavy-quark limit Ri{w) = R2{w) — 1 independently of w. 

To extract the functions hAi{w)^ Ri{w) and R2{w) from experimental data 
is a complicated task. However, HQET-based calculations suggest that the w 
dependence of the fopxi-factor ratios, which is induced by symmetry-breaking 
effects, is rather mildE3. Moreover, the form factor hAiiw) is expected to have 
a nearly linear shape over the accessible w range. This motivates to introduce 
three parameters g^-^^, i?i and R2 by 



hAi{w) « 
Ri{w) « Ri , 



1- g\^{w-l) 
R2{w) « i?2 



(65) 



where = 0.91 ± 0.03 from (||). The CLEO Collaboration has extracted 
these three parameters from an analysis of the angular distributions in _B ^ 
D*iP decays l3. The results are 



0.91 ±0.16, i?i = 1.18 ±0.32, 



R2 



0.71 ±0.23. 



(66) 



Using the HQET, one obtains an essentially model-independent prediction for 
the symmetry-breaking corrections to whereas the corrections to i?2 are 
somewhat model dependent. To good approximation 



Ri 



i?2 



1 



A 



1 - K- 



A 

2'mr 



2mc 
0.8 ±0.2 



1.3 ±0.1, 



(67) 
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with K « 1 from QCD sum rule^ 



Jlere A is the "binding energy" as define 
1 as weh as phenomenological analysesHi 



in (^. Theoretical calculationsBifi 
suggest that A « 0.45-0.65 GeV is the appropriate value to be used in one-loop 
calculations. A quark- model calculation of i?i and i?2 gives results similar to 
the HQET predictions l3: Ri ~ 1.15 and R2 « 0.91. The experimental data 
confirm the theoretical prediction that i?i > 1 and i?2 < 1, although the errors 
are still large. 

Heavy-quark symmetry has also been tested by comparing the form factor 
^{w) in i? — > D*£D decays with the corresporiding form factor G{w) governing 
B D £D decays. The theoretical predictionUo 

= 1.08 ± 0.06 (68) 

compares well with the expeiipiental results for this ratio: 0.99 ±0.19 reported 
by the ClfEO Collaboration^ and 0.87±0.30 reported by the ALEPH CoUab- 
orationEEI. In these analyses, it has also been tested that within experimental 
errors the shape of the two form factors agrees over the entire range of w values. 

The results of the analyses described above are very encouraging. Within 
errors, the experiments confirm the HQET predictions, starting to test them 
at the level of symmetry-breaking corrections. 



4 Inclusive Decay Rates 

Inclusive decay rates determine the probability of the decay of a particle into 
the sum of all possible final states with a given set of global quantum num- 
bers. An example is provided by the inclusive semi-leptonic decay rate of the 
B meson, T{B — > X £i>), where the final state consists of a lepton-neutrino pair 
accompanied by any number of hadrons. Here we shall discuss the theomti 
cal description of inclusive decays of hadrons containing a heavy quark 
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From a theoretical point of view such decays have two advantages: first, bound- 
state effects related to thainitial state, such as the "Fermi motion" of the heavy 
quark inside the hadronE23ll23, can be accounted for in a systematic way us- 
ing the heavy-quark expansion; secondly, the fact that the final state consists 
of a sum over many hadronic channels eliminates bound-state effects related 
to the properties of individual hadrons. This second feature is based on the 
hypothesis of quark-hadron duality, which is an important concept in QCD 
phenomenology. The assumption of duality is that cross sections and decay 
rates, which are defined in the physical region (i.e. the region of time-like mo- 
menta), are calculable in QCD after a "smearing" or "averaging" procedure 
has been applied t23. In semi-leptonic decays, it is the integration over the 
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lepton and neutrino phase space that provides a smearing over the invariant 
hadronic mass of the final state (so-called global duality) . For non-leptonic de- 
cays, on the other hand, the total hadronic mass is fixed, and it is only the fact 
that one sums over many hadronic states that provides an averaging (so-called 
local duality) . Clearly, local duality is a stronger assumption than global dual- 
ity. It is important to stress that quark-hadron duality cannot yet be derived 
from first principles; still, it is a necessary assumption for many applications 
of QCD. The validity of glnbal duality has been tested experimentally using 
data on hadronic r decaysE23. 

Using the optical theorem, the inclusive decay width of a hadron Hf, con- 
taining a b quark can be written in the form 

T{Hb^X) = -^lm{Hb\T\Hb), (69) 

where the transition operator T is given by 

T^ijd''xT{C,s{x),C,s{m- (70) 

Inserting a complete set of states inside the time-ordered product, we recover 
the standard expression 

T{Hb ^ X) ^ -^Y.(2n)H\pH - Px)\{X\C,s\Hb)? (71) 

for the decay rate. For the case of semi-leptonic and non-leptonic decays, Ccs 
is the effective weak LagimLgiaii ^iven in (^, which in practice is corrected 



for short-distance effects Sii3'tila~Liill arising from the exchange of gluons with 



virtualities between mw and m?,. If some quantum numbers of the final states 
X are specified, the sum over intermediate states is to be restricted appropri- 
ately. In the case of the inclusive semi-leptonic decay rate, for instance, the 
sum would include only those states X containing a lepton-neutrino pair. 

In perturbation theory, some contributions to the transition operator are 
given by the two-loop diagrams shown on the left-hand side in Fig. |[ Because 
of the large mass of the b quark, the momenta flowing through the internal 
propagator lines are large. It is thus possible to construct an OPE for the 
transition operator, in which T is represented as a series of local operators 
containing the heavy-quark fields. The operator with the lowest dimension, 
d = 3, is bh. It arises by contracting the internal lines of the first diagram. 
The only gauge-invariant operator with dimension 4 is bipb; however, the 
equations of motion imply that between physical states this operator can be 
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Figure 9: Perturbative contributions to the transition operator T (left), and the correspond- 
ing operators in the OPE (right). The open squares represent a four-fermion interaction of 
the effective Lagrangian £cff i while the black circles represent local operators in the OPE. 

replaced by mi,bb. The first operator that is different from bb has dimension 5 
and contains the gluon field. It is given by b ggCFf^iyG^'^b. This operator arises 
from diagrams in which a gluon is emitted from one of the internal lines, such 
as the second diagram shown in Fig. ^. For dimensional reasons, the matrix 
elements of such higher-dimensional operators are suppressed by inverse powers 
of the heavy-auaxk mass. Thus, any inclusive decay rate of a hadron Hh can 
be written asEZl^ 

^^^^^^^^-^M^Y^^^^)" + '^ + ---|^ (72) 

where the prefactor arises naturally from the loop integrations, are calcu- 
lable coefficient functions (which also contain the relevant CKM matrix ele- 
ments) depending on the quantum numbers / of the final state, and {0)h are 
the (normalized) forward matrix elements of local operators, for which we use 
the short-hand notation 

{0)h = m O \H,) . (73) 

In the next step, these matrix elements are systematically esBaadcd in 
powers of 1/mb, using the technology of the HQET. The result isEa^EZlej 

{bb), = l- ''"^''''\'f^^''''K 0il/ml), 

{bg,a^,G'"'b)H = 2^il{Hi,) + 0(l/mb) , (74) 
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where we have defined the HQET matrix elements 

^il{H,) = {Ht{v)\ K (iDf K \Hb{v)) , 

MG(^b) = TT^ {Hk{v)\ bv^<J^.G'^''b. \Hb{v)) . (75) 

Here (iD)^ = {iv ■ D)^ -— (iD)^; in the rest frame, this is the square of the 
operator for the spatial momentum of the heavy quark. Inserting these results 
into ( l?^ ) yields 

(76) 

It is instructive to understand the appearance of the "kinetic energy" contri- 
bution /x^, which is the gauge-covariant extension of the square of the 6-quark 
momentum inside the heavy hadron. This contribution is the field-theory ana- 
logue of the Lorentz factor (1 — v^Y^^ ~ 1 — fc^/2m^, in accordance with the 
fact that the lifetime, r = l/F, for a moving particle increases due to time 
dilation. 

The main result of the heavy-quark expansion for inclusive decay rates 
is the observation that the free quark decay (i.e_the parton model) provides 
the first term in a systematic l/rrib expansion tj. For dimensional reasons, 
the corresponding rate is proportional to the fifth power of the 6-quark mass. 
The non-perturbative corrections, which arise from bound-state effects inside 
the B meson, are suppressed by at least two powers of the heavy-quark mass, 
i.e. they are of relative order {Aqcu / rrib)'^ . Note that the absence of first- 
order power corrections is a consequence of the equations of motion, as there 
is no independent gauge-invariant operator of dimension 4 that could appear 
in the OPE. The fact that bound-state effects in inclusive decays are strongly 
suppressed exnlai.ys a posteriori the success of the parton model in describing 
such processes tilHli^. 

The hadronic matrix elements appearing in the heavy-quark expansion 
( [76| ) can be determined to some extent from the known masses of heavy hadron 
states. For the B meson, one finds that 

Ai^(S) = -Ai = (0.3 ± 0.2) GeV^ , 

filiB) = 3A2 w 0.36 GeV^ , (77) 



where Ai and A2 are the parameters appearing in the mass formula (|32). For 
the ground-state baryon A^, in which the light constituents have total spin 
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zero, it follows that 



Mg(A6) = 0, 



while the matrix element /i^(Ab) obeys the relation 



{mA^ - rriAj - {ruB - mo) ^ fil{B) - ^l{Ab) 



1 



2mr 



2mf, 



(78) 

+ 0(l/m|), 
(79) 

where ttib and m^i denote the spin-averaged masses introduced in connection 
with The above relation implies 

filiB) - filiAb) = (0.01 ± 0.03) GeV' . (80) 

What remains to be calculated, then, is the coefhcient functions for a given 
inclusive decay channel. 

To illustrate this general formalism, we discuss as an example the de- 
termination of \Vcb\ from inclusive semi-leptonic B decays. In this xiasa^the 
short-distance coefficients in the general expression (|76|) are given bytZlLJ 



'6\Vcb\Hl 



8x^ + 8x^ 
xY ■ 



12x^ \nx^ +0{as 



(81) 



Here x = nic/mb, and rrib and nic are the masses of the b and c quarks, 
defined to a gi ven order in perturbation theoryEj. The 0{as) terms in i 



,SL 



known exactly 



U3, and reliable estimates exist for the 0{a^) corrections Eif 
The theoretical uncertainties in this determination of \Vcb\ are quite different 
from those entering the analysis of exclusive decays. The main sources are the 
dependence on the heavy-quark masses, higher-order perturbative corrections, 
and above all the assumption of global quark-hadron duality. A conservatisifi 
estimate of the total theoretical error on the extracted value of \Vcb\ yields 
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\Vcb\ = (0.040±0.003) 



BsL 


1/2 


1.6 ps 


1/2 


10.5% 









(40±lcxp±3th)xl0-3. (82) 



The value of |Vcf,| extracted from the inclusive semi-leptonic width is in ex- 
cellent agreement with the value in (^3|) obtained from the analysis of the 
exclusive decay B — > D*£ v. This agreement is gratifying given the differences 
of the methods used, and it provides an indirect test of global quark-hadron 
duality. Combining the two measurements gives the final result 

iVcbl = 0.039 ± 0.002. (83) 

After Vud and Vus, this is the third-best known entry in the CKM matrix. 
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5 Rare B Decays and Determination of the Weak Phase 7 



The main objectives of the B factories are to explore the physics of CP vio- 
lation, to determine the flavor parameters of the electroweak theory, and to 
probe for physics beyond the Standard Model. This will test the CKM mecha- 
nism, which predicts that all CP violation results from a single complex phase 
in the quark mixing matrix. Facing the announcement of evidence fiju^ CP 
asymmetry in the decays B — > J/ipRs by the CDF Collaboration ti3, the 
of direct CP violation in X —> tttt decays by the KTeV and NA48 
and the successful start of the B factories at SLAC and KEK, 
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groups t' 



the year 1999 has been an important step towards achieving this goal. 




(0,0) 



(1,0) 



Figure 10: The rescaled unitarity triangle representing the relation 1- 



: 0. 



The apex is determined by the Wolfenstein parameters (p, r;). The area of the triangle is 
proportional to the strength of CP violation in the Standard Model. 



The determination of the sides and angles of the "unitarity triangle" 
V*^Vud + V*ijVcd+VflVtd = depicted in Fig. |l^ plays a central role in the B fac- 
tory program. Adopting the standard phase conventions for the CKM matrix, 
only the two smallest elements in this relation, V*ij and Vtd, have non- vanishing 
imaginary parts (to an excellent approximation). In the Standard Model the 
angle (3 =^ — arg(Vfd) can be determined in a theoretically clean way by mea- 
suring the mixing-induced CP ^asymmetry in the decays B — > J/ipKs- The 



preliminary CDF result imphcsEia sin 2/3 = 0.79ljj:44. The angle 7 = arg(y„*J, 
or pniiiya.lently the combination a = 180° — /5 — 7, is much harder to deter- 
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Recently, there has been significant progress in the theoretical under- 
standing of the hadronic decays B ttK, and methods have been developed 
to extract information on 7 from rate measurements for these processes. Here 
we discuss the charged modes B^ ttK, which from a theoretical perspective 
are particularly clean. 
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In the Standard Model, the main contributions to the decay amphtudes 
for the rare processes B ttK are due to the penguin-induced flavor-changing 
neutral current (FCNC) transitions b — > sqq, which exceed a small, Cabibbo- 
suppressed h — > uus contribution from VF-boson exchange. The weak phase 7 
enters through the interference of these two ( "penguin" and "tree" ) contribu- 
tions. Because of a fortunate interplay of isospin, Fierz and flavor symmetries, 
the theoretical description of the charged modes — > -n-K \r vprv clean de- 



spite the fact that these are exclusive non-leptonic decaysEijil^. Without any 
dynamical assumption, the hadronic uncertainties in the description of the in- 
terference terms relevant to the determination of 7 are of relative magnitude 
O(A^) or 0{esu(3)/-^c), where A = sinOc ~ 0.22 is a measure of Cabibbo sup- 
pression, esu(3) ~ 20% is the typical size of SU(3) breaking, and the factor l/iVc 
indicates that the corresponding terms vanish in the factorization approxima- 
tion. Factorizable SU(3) breaking can be accounted for in a straightforward 
way. 

Recently, the accuracy of this description has been further improved when 
it was shown that non-leptonic B decays into two light mesons, pSuch as i3 ^ 
ttK and B tttt, admit a systematic heavy-quark expansion EHa. To lead- 
ing order in 1/mf,, but to all orders in perturbation theory, the decay am- 
plitudes for these processes can be calculated from first principles without 
recourse to phenomenological models. The QCD factorization theorem proved 
in Ref. 122 iuaproves upon the phenomenological approach of "generalized fac- 
torization" which emerges as the leading term in the heavy-quark limit. 
With the help of this theorem, the irreducible theoretical uncertainties in the 
description of the B^ ttK decay amplitudes can be reduced by an extra 
factor of 0(l/m{,), rendering their analysis essentially model independent. As 
a consequence of this fact, and because they are dominated by FCNC tran- 
sitions, the decays R'*' r^^^ oS^r a sensitive probe to physics beyond the 



Standard Model 
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much in the same way as the "classical" FCNC 



processes B — > ^^7 or B ^ t^l 



5.1 Theory of B^ ttK Decays 

The hadroni c de cays B — s- ttK are mediated by a low-energy effective weak 
Hamiltonian E£3, whose operators allow for three different classes of flavor 
topologies: QCD penguins, trees, and electroweak penguins. In the Standard 
Model the weak couplings associated with these topologies are known. From 
the measured branching ratioS|-oiie can deduce that QCD penguins dominate 
the B ttK decay amplitudesllHa, wkereas trees and electroweak penguins are 
subleading and of a similar strengthll23. The theoretical description of the two 
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charged modes tt^K^ and tt'^K^ exploits the fact that the am- 

phtudes for these processes differ in a pure isospin amphtude, ^3/2, defined as 
the matrix element of the isovector part of the effective Hamiltonian between a 
B meson and the ttK isospin eigenstate with I = ^- In the Standard Model the 
parameters of this amplitude are determined, up to an overall strong phase (f), 



in the limit of SU(3) flavor symmetry till. Using the QCD factorization theorem 



the SU(3)-breaking corrections can be calculated in a model-independent way 
up to non-factorizable terms that are power-suppressed in 1 /mb and vanish in 
the heavy-quark limit. 

A convenient parameterization of the non-leptoniC|depay amplitudes ^+0 
AiB+ TT+K") and Ao+ = -V2A{B+ n"K+) isli 



^+0 = P(l-£ae^^e"'), 



Ao+ = P 



l-e„e^^e"'-£3/2e'^e'^-<5Ew) , (84) 



where P is the dominant penguin amplitude defined as the sum of all terms in 
the B^ tt^K'^ amplitude not proportional to e*''', 77 and are strong phases, 
and Eq, £3/2 and few are real hadronic parameters. The weak phase 7 changes 
sign under a CP transformation, whereas all other parameters stay invariant. 

Based on a naive quark-diagram analysis one would not expect the B^ — > 
TT'^K'^ amplitude to receive a contribution from b uus tree topologies; how- 
ever, such a contribution xanp-be induced through final-state rescattering or 



annihilation contributionslH2Hi2£l. They are parameterized by Sa = 0(A ). In 



the heavy-quark limit this parameter can be calculated and is found to be very 



small IHJ: Sa ~ —2%. In the future, it will be possible to put upper and lower 
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bounds on Sa by comparing the CP-averaged branching ratios for the decays 
B^ TT^K^ and B^ K^K^. Below we assume \ea\ < 0.1; however, our 
results will be almost insensitive to this assumption. 

The terms proportional to £3/2 in ( p^ parameterize the isospin amplitude 
A3/2. The weak phase e^'^ enters through the tree process b — > uus, whereas 
the quantity few describes the effects of electroweak penguins. The parame- 
ter £3/2 measures the relative strength of tree and QCD penguin contributions. 
Information about it can be derived by using SU(3) flavor symmetry to relate 
the tree contribution to the isospin amplitude to the corresponding con- 
tribution in the decay 5+ tt+tt^. Since the final state Tr+Tr*^ has isospin 
1 = 2, the amplitude for this process does not receive any contribution from 
QCD penguins. Moreover, electroweak penguins in b —^ dqq transitions are 
negligibly small. We define a related parameter £3/2 by writing 



'3/2 



£3/2 a/1 - 2£a COS ?7 cos 7 -I- el , (85) 
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so that the two quantities agree in the Umit Ea 0. In the SU(3)_liaiit this 
new parameter can be determined experimentaUy form the relationlHj 



-3/2 



Ri 



2B(B±^^^r±7rO) 



1/2 



(86) 



SU(3)-breaking corrections are described by the factor i?i = 1.22 ±0.05, which 
can be calculated in a model-indfipendent way using the QCD factorization 
theorem for non-leptonic decaysEHI. The quoted error is an estimate of the the- 
oretical uncertainty due to correctit^ns of 0(-^^). Using prehminary data re- 



ported by the CLEO Collaboration 
branching ratios in (B6h, we obtain 
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to evaluate the ratio of the CP-averaged 



^3/2 



0.21 ±0.06cxp± 0.01th. 



(87) 



With a better measurement of the branching ratios the uncertainty in 63/2 will 
be reduced significantly. 
Finally, the parameter 



3EW 



i?2 



v*v 

^cb^cs 



ub 



= (0.64 ±0.09) 



87r siv? 9w 
0.085 

\VublVc 



1 



3 In Xt 



cb\ 



with Xt — {mt/mw)'^, describes the ratio of electroweak penguin and tree 
contributions to the isospin amplitude A^if the SU(3) limit it is calculable 

.ll£,L3s| 



in terms of Standard Model parameters 



SU(3)-breaking as wpII 



s small 

electromagnetic corrections are accounted for by the quantity F^^H^^I R2 = 
0.92 ±0.09. The error quoted in ( |8^ ) includes the uncertainty in the top-quark 
mass. 

Important observables in the study of the weak phase 7 are the ratio of 
the CP-averaged branching ratios in the two ttK decay modes. 



i?* — 



^ 0.75 ±0.28, 



2B(S± ^ ttOR^) 
and a particular combination of the direct CP asymmetries 



A = 



0.52 ±0.42. 



(89) 



(90) 



The experimental values of these quantiiies are derived using preliminary data 



reported by the CLEO Collaboration liiia. The theoretical expressions for i?* 
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and A obtained using the parameterization in ( p^ are 

R^^ = 1 + 263/2 cos (/-(Jew -cos 7) 

+ ^3/2(1 - 2(5ew C0S7 + (^|w) + 0(e3/2ea) , 

A = 2£3/2sin7sin0 + O(£3/2£Q). (91) 

Note that the rescattering effects described by Sa are suppressed by a factor 
of £3/2 and thus reduced to the percent level. Explicit expressions for these 
contributions can be found in Ref. 121. 



5.2 Lower Bound on 7 and Constraint in the {p,ri) Plane 

There are several strategies for exploiting the above relations. From a mea- 
surement of the ratio i?, alone a bound on cos 7 can be derived, implying a 
non-trivial constraint on jthe Wolfenstein parameters p and fj defining the apex 
of the unitarity triangle 1113. Only CP-averaged branching ratios are needed 
for this purpose. Varying the strong phases (j) and rj independently we first 
obtain an_upper bound on the inverse of i?* . Keeping terms of linear order in 
£a yields 
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i?* ^ < (1 +£3/2 Ifew - COS7I) -|-£3/2sin^7-f 2£3/2|£Q|sin^7. (92) 

Provided i?, is significantly smaller than 1, this bound implies an exclusion 
region for cos 7 which becomes larger the smaller the values of /L_and £3/2 
are. It is convenient to consider instead of i?, the related quantity^ 



Xr^^-^ = 0.72 ±0.98cxp± 0.03th. (93) 

£3/2 



Because of the theoretical factor i?i entering the definition of £3/2 in (ISq ) this 
is, strictly speaking, not an observable. However, the theoretical uncertainty 
in Xr is so much smaller than the present experimental error that it can be 
ignored for all practical purposes. The advantage of presenting our results in 
terms of Xr rather than is that the leading dependence on £3/2 cancels 
out, leading to the simple bound \Xr\ < \Sew — COS7I -I- 0(£3/2, £a). 

In Fig. ^ we show the upper bound on Xr as a function of I7I, obtained 
by varying the input parameters in the intervals 0.15 < £3/2 < 0.27 and 0.49 < 
Sew < 0.79 (corresponding to using \Vub/Vcb\ = 0.085 ± 0.015 in (H)). Note 
that the effect of the rescattering contribution parameterized by Sa is very 
small. The gray band shows the current value of Xr, which still has too large 
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Figure 11: Theoretical upper bound on the ratio Xr versus I7I for ea = 0.1 (soHd Hne) 
and ea = (dashed Hne). The horizontal line and band show the current experimental value 

with its 1(T variation. 



an error to provide any useful information on 7. The situation may change, 
however, once a more precise measurement of Xn will become available. For 
instance, if the current central value = 0.72 were confirmed, it would imply 
the bound I7I > 75°, marking a significant improvement over the indirect limit 
I7I > 37° inferred from the globaJ analysis of the unitarity triangle including 



So far, we have used the inequality ( JQ^ ) to derive a lower bound on |7|. 
However, a large part of the uncertainty in the value of (Jewi and thus in the 
resulting bound on I7I, comes from the present large error on \Vub\- Since this 
is not a hadronic uncertainty, it is appropriate to separate it and turn into 
a constraint on the Wolfenstein parameters p and fj. To this end, we use that 
COS7 = pI \/ + Tf by definition, and few = (0.24 ± 0.03)/\/p^ + ff- from 
(^). The solid lines in Fig. |l^ show the resulting constraint in the (p, 77) plane 
obtained for the representative values Xji = 0.5, 0.75, 1.0, 1.25 (from right to 
left), which for £3/2 = 0.21 would correspond to i?, = 0.82, 0.75, 0.68, 0.63, 
respectively. Values to the right of these lines are excluded. For comparison, 
the dashed circles show the constraint arising from the measurement of the ratio 
\Vub/Vcb\ = 0.085 ± 0.015 in semi-leptonic B decays, and the dashed-dotted 
line shows the bound implied b y th e present experimental limit on the mass 
difference Arn^ in the Bg systemEij. Values to the left of this line are excluded. 
It is evident from the figure that the bound resulting from a measurement of 
the ratio Xn in B^ ttK decays may be very non-trivial and, in particular, 
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Figure 12: Theoretical constraints on the Wolfenstein parameters (p, f?) impUed by a mea- 
surement of the ratio X/j in — > nK decays (sohd hnos), semi-leptonic B decays (dashed 
circles), and B^^-B^ mixing (dashed-dotted line). 



may eliminate the possibility that 7 = 6. The combination of this bound 
with information from semi-leptonic decays and B~B mixing alone would then 
determine the Wolfenstein parameters p and fj within narrow ranges|| and in 
the context of the CKM model would prove the existence of direct CP violation 
in B decays. If one is more optimistic, one may even hope that in the future 
the constraint from B — > nK decays may become incompatible with the bound 
from Bg-Bs mixing, thus indicating New Physics beyond the Standard Modelfl 



5.3 Extraction of j 

Ultimately, the goal is of course not only to derive a bound on 7 but to de- 
termine this parameter directly from the data. This requires to fix the strong 
phase in which can be achieved either through the measurement of a 
CP asymmetry or with the help of theory. A strategy for an experimental 
determination of 7 from B^ nK decays has been suggested ixi_Ref. 120. It 
generalizes a method proposed by Gronau, Rosner and LondonlllJ to include 
the effects of electroweak penguins. The approach has later been refined to^ 
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count for rescattering contributions to the B^ — > tt^K^ decay amplitudes 
Before discussing this method, we will first illustrate an easier strategy for a 
theory-guided determination of 7 based on the QCD factorization theorem for 



'^An observation of CP violation, such as the measurement of in K-K mixing or sin 2/3 
in B — > J/ip Kg decays, is however needed to fix the sign of fj. 

''At the time of writing, the bound from Bs-Bs mixing is being pushed further to the 
right, making such a scenario a tantalizing possibility. 
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non-lcptonic decaysll£3. This method does not require any measurement of a 
CP asymmetry. 

Theory- guided determination: 

In the previous section the theoretical predictions for the non-leptonic B — s- tiK 
decay amphtudes obtained using the QCD factorization theorem were used in 
a minimal way, i.e. only to calculate the size of the SU(3)-breaking effects 
parameterized by i?i and i?2 in ( p6| ) and (|8|). The resulting bound on 7 and 
the corresponding constraint in the (p, fj) plane are therefore theoretically very 
clean. However, they are only useful if the value of X]^ is found to be larger 
than about 0.5 (see Fig. ^l|), in which case values of I7I below 65° are excluded. 
If it would turn out that Xji < 0.5, then it is possible to satisfy the inequality 
( p^ ) also for small values of 7, however, at the price of having a very large 
strong phase, cj) « 180°. But this aiassibility can be discarded based on the 
model-independent prediction thatli23 

= 0[as{mb), AQCD/mt] . (94) 

A direct calculation of this phase to leading power in l/m^ yieldslH3 cj) w —11°. 
Using the fact that (f> is parametrically small, we can exploit a measurement 
of the ratio Xji to obtain a determination of I7I - corresponding to an allowed 
region in the (p, fj) plane - rather than just a bound. This determination is 
unique up to a sign. Note that for small values of (p the impact of the strong 
phase in the expression for in ( pl| ) is a second-order effect. As long as 
|(/)| ^ y^2Aiy^J¥y^, the uncertainty in coscj) has a much smaller effect than 
the uncertainty in £3/2. With the present value of £3/2 this is the case as long 
as 1^1 ^ 43°. We believe it is a safe assumption to take < 25° (i.e. more 
than twice the value obtained to leading order in 1/mf,), so that coscj) > 0.9. 

Solving the equation for _R» in ( pl[ ) for cos 7, and including the corrections 
of 0{ea), we find 

Xji + ie3/2(X2 - 1 -t- Sly^) Ea cos77sin27 
cos7 = (5ew —— — 7 + — — — - — , (95) 

COS0 -I- £3/20EW COS0 -I- £3/20EW 

where we have set coscj) = 1 in the numerator of the 0(£a) term. Using the 
QCDpiaptorization theorem one finds that £a cos 77 ~ —0.02 in the heavy-quark 
limit HHI, and we assign a 100% uncertainty to this estimate. In evaluating the 
result ( |95| ) we scan the parameters in the ranges 0.15 < £3/2 < 0.27, 0.55 < 
few < 0.73, -25° < </) < 25°, and -0.04 < £aCos?7sin27 < 0. Figure |l| shows 
the allowed regions in the (p, ff) plane for the representative values Xfi = 0.25, 
0.75, and 1.25 (from right to left). We stress that with this method a useful 
constraint on the Wolfenstein parameters is obtained for any value of Xr. 
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Figure 13: Allowed regions in the (p, r;) plane for fixed values of Xn, obtained by varying 
all theoretical parameters inside their respective ranges of uncertainty, as specified in the 
text. The sign of f; is not determined. 



Model-independent determination: 

It is important that, once more precise data on nK decays will become 

available, it will be possible to test the prediction of a small strong phase 
(j) experimentally. To this end, one must determine the CP asymmetry A 
defined in ( pO| ) in addition to the ratio R^. From it follows that for fixed 
values of £3/2 and Sew the quantities R^: and A define contours in the (7, 0) 
plane, whoaaJxitfirsections determine the two phases up to possible discrete 



ambiguities lH3'E£ll. Figure ^ shows these contours for some representative 
values, assuming £3/2 = 0.21, Sew = 0.64, and Sa = 0. In practice, including 
the uncertainties in the values of these parameters changes the contour lines 
into contour bands. Typicallv^ 



determination of 7 of about 
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tjie spread of the bands induces an error in the 
10°. In the most general case there are up 
to eight discrete solutions for the two phases, four of which are related to the 
other four by a sign change (7, 0) (—7,-0). However, for typical values of 
i?* it turns out that often only four solutions exist, two of which are related 
to the other two by a sign change. The theoretical prediction that cj) is small 
implies that solutions should exist where the contours intersect close to the 
lower portion in the plot. Other solutions with large are strongly disfavored. 
Note that according to (^l|) the sign of the CP asymmetry A fixes the relative 
sign between the [twp phases 7 and ((>. If we trust the theoretical prediction 



that (j) is negativeliiil, it follows that in most cases there remains only a unique 



solution for 7, i.e. the CP-violating phase 7 can be determined without any 
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Figure 14: Contours of constant ij, ("hyperbolas") and constant |j4| ("circles") in the 
(I7I, \<j>\) plane. The sign of the asymmetry A determines the sign of the product sin7sin(/). 
The contours for ij* refer to values from 0.6 to 1.0 in steps of 0.1, those for the asymmetry 
correspond to 5%, 15%, and 25%, as indicated. 



discrete ambiguity. 

Consider, as an example, the hypothetical case where i?* = 0.8 and A = 
-15%. Figure |l| then allows the four solutions where (7,(/i) w (±82°,=f21°) 
or (±158°, =f78°). The second pair of solutions is strongly disfavored because 
of the large values of the strong phase 0. From the first pair of solutions, the 
one with (j) « —21° is closest to our theoretical expectation that (j) ~ 
hence leaving 7 « 82° as the unique solution. 



6 Sensitivity to New Physics 

In the presence of New Physics the theoretical description of — > nK decays 
becomes more complicated. In particular, new CP-violating contributions to 
the decay amplitudes may be induced. A detailed analysis of such effects has 



been presented inlirH. A convenient and completely general parameterization 



of the two amplitudes in ( p4[ ) is obtained by replacing 

P^P' , Ea e'^e*" ^ ip e"^" , few ^ a e'"^" + ih e"^" , (96) 
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where p, a, b are real hadronic parameters, and (jjp, (pa, (f'b are strong phases. 
The terms ip and ib change sign under a CP transformation. New Physics 
effects parameterized by P' and p are isospin conserving, while those described 
by a and b violate isospin symmetry. Note that the parameter P' cancels in all 
ratios of branching ratios and thus does not affect the quantities i?* and Xn 
as well as any CP asymmetry. Because the ratio i?* in ( |8^) would be 1 in the 
limit of isospin symmetry, it is particularly sensitive to isospin-violating New 
Physics contributions. 

New Physics can affect the bound on 7 derived from (|9^) as well as the 
extraction of 7 using the strategies discussed above. We will discuss these two 
possibilities in turn. 



6.1 Effects on the Bound on 7 

The upper bound on R^^ in ( p^ and the corresponding bound on Xn shown 
in Fig. |l^ are model-independent results valid in the Standard Model. Note 
that the extremal value of is such that \Xji\ < (1 + Sew) irrespective 
of 7. A. va]ue.jaL|Xji;| exceeding this bound would be a clear signal for New 
Physics llBllHilEZl. 

Consider first the case where New Physics may induce arbitrary CP- 
violating contributions to the B — ^ t:K decay amplitudes, while preserving 
isospin symmetry. Then the only change with respect to the Standard Model 
is that the parameter p may no longer be as small as 0(£a)- Varying the 
strong phases (p and (j)p independently, and allowing forj-an arbitrarily large 
New Physics contribution to p, one can derive the boundEHU 

\Xb.\ < \/i - 2(5ew C0S7 + < 1 + Sew ■ (97) 

The extremal value is the same as in the Standard Model, i.e. isospin-conserving 
New Physics effects cannot lead to a value of \Xji\ exceeding {1 + Sew)- For in- 
termediate values of 7 the Standard Model bound on Xj^ is weakened; but even 
for large p = 0(f), corresponding to a significant New Physics contribution to 
the decay amplitudes, the effect is small. 

If both isospin-violating and isospin-conserving New Physics contributions 
are present and involve new CP-violating phases, the analysis becomes more 
complicated. Still, it is possible to derive model- independent bounds-on Xn. 
Allowing for arbitrary values of p and all strong phases, one obtainsll23 

\Xr\ < V(|a| + |cos7|)2 + (|6| + |sin7|)2 



< f + Va^ + &2 < _^Xr, (98) 

£3/2 
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where the last inequahty is relevant only in cases where \J a? ^ 1. The 
important point to note is that with isospin-violating New Physics contribu- 
tions the value of \Xji\ can exceed the upper bound in the Standard Model by 
a potentially large amount. For instance, if is twice as large as in the 

Standard Model, corresponding to a New Physics contribution to the decay 
amplitudes of only 10-15%, then \X^\ could be as large as 2.6 as compared 
with the maximal value 1.8 allowed (for arbitrary 7) in the Standard Model. 
Also, in the most general case where h and p are non-zero, the maximal value 
jXfll can take is no longer restricted to occur at the endpoints 7 = 0? or 180°, 
which are disfavored by the global analysis of the unitarity triangleEij. Rather, 
\Xy{\ would take its maximal value if |tan7| = \p\ = |fc/a|. 

The present experimental value of in (|9^) has too large an error to 
determine whether there is any deviation from the Standard Model. If Xr 
turns out to be larger than 1 (i.e. at least one third of a standard deviation 
above its current central value), then an interpretation of this result in the 
Standard Model would require a large value I7I > 91° (see Fig. |ll|), which 
would be difficult to accommodate in view of the upper bound implied by the 
experimental constraint on Bg-Bg mixing, thus providing evidence for New 
Physics. If Xji > 1.3, one could go a stefL,further and conclude that the New 
Physics must necessarily violate isospin 
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6.2 Effects on the Determination of j 

A value of the observable i?* violating the bound ( |9^ ) would be an exciting 
hint for New Physics. However, even if a future precise measurement will give 
a value that is consistent with the Standard Model bound, B^ nK decays 
provide an excellent testing ground for physics beyond the Standard Model. 
This is so because New Physics may cause a significant shift in the value of 7 
extracted using the strategies discussed earlier, leading to inconsistencies when 
this value is compared with other determinations of 7. 

A global fit of the unitarity triangle combining information from semi- 
leptonic B decays, B~B mixing, CP violation in the kaon system, and mixing- 
induced CP violation in B —^ J ftp Ks decays provides information qti 7 which 



in a few years will determine its value within a rather narrow range tij. Such 
an indirect determination could be complemented by direct measurements of 7 
using, e.g., B DK^*^ decays, or using the triangle relation 7 — 180° — a — (3 
combined with a measurement of a. We will assume that a discrepancy of 
more than 25° between the "true" 7 = arg(yj*j) and the value ^t^k extracted 
in B^ ttK decays will be observable after a few years of operation at the B 
factories. This sets the benchmark for sensitivity to New Physics effects. 
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Figure 15: Contours of constant Xji versus 7 and the parameter a, assuming 7 > 0. The 
horizontal band shows the value of a in the Standard Model. 



In order to illustrate how big an effect New Physics could have on the 
extracted value of 7, we consider the simplest case where there are no new 
CP-violating couplings. Then all New Physics contributions in ( |9^ ) are param- 
eterized by the single parameter cnp = a — Sew- A more general discussion 
can be found in Ref. 127. We also assume for simplicity that the strong phase 
(j) is small, as suggested by (|4|). In this case the difference between the value 
Jttk extracted from ttK decays and the "true" value of 7 is to a good 

approximation given by 

cos 77r/<- — cos 7 — cnp • (99) 

In Fig. |l^ we show contours of constant versus 7 and a, assuming without 
loss of generality that 7 > 0. Obviously, even a moderate New Physics contri- 
bution to the parameter a can induce a large shift in 7. Note that the present 
central value of Xr « 0.7 is such that values of a less than the Standard Model 
result a ~ 0.64 are disfavored, since they would require values o£-q^-|exceeding 
100°, in conflict with the global analysis of the unitarity triangle E 
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6.3 Survey of New Physics models 

In Ref. 127, we have explored how New Physics could affect purely hadronic 
FCNC transitions of the type b — > sqq focusing, in particular, on isospin vio- 
lation. Unlike in the Standard Model, where isospin-violating effects in these 



44 



Table 1: Maximal contributions to a^p and shifts in 7 in extensions of the Standard Model. 
For the case of supcrsymmctric models with R-parity the first (second) row corresponds 
to maximal right-handed (left-handed) strange-bottom squark mixing. For the two-Higgs- 
doublet models we take mjj+ > 100 GeV and tan /3 > 1. 



Model 




FCNC Z exchange 

extra Z' boson 

SUSY without R-parity 


2.0 180° 
14 180° 
14 180° 


SUSY with R-parity 


0.4 25° 
1.3 180° 


two-Higgs-doublet models 
anomalous gauge-boson couplings 


0.15 10° 
0.3 20° 



processes are suppressed by electroweak gauge couplings or small CKM matrix 
elements, in many New Physics scenarios these effects are not parametrically 
suppressed relative to isospin-conserving FCNC processes. In the language of 
effective weak Hamiltonians this implies that the Wilson coefficients of QCD 
and electroweak penguin operators are of a similar magnitude. For a large 
class of New Physics models we found that the coefficients of the electroweak 
penguin operators are, in fact, due to "trojan" penguins, which are neither 
related to penguin diagrams nor of electroweak origin. 

Specifically, we have considered: (a) models with tree-level FCNC cou- 
plings of the Z boson, extended gauge models with an extra Z' boson, su- 
persymmetric models with broken R-parity; (b) supersymmetric models with 
R-parity conservation; (c) two-Higgs-doublet models, and models with anoma- 
lous gauge-boson couplings. Some of these models have also been investigated 
in Refs. 125 and 126. In case (a), the electroweak penguin coefficients can be 
much larger than in the Standard Model because they are due to tree-level pro- 
cesses. In case (b), these coefficients can compete with the ones of the Standard 
Model because they arise from strong- interaction box diagrams, which scale 
relative to the Standard Model like (as/Q;)(TO^y/m|ugY). In models (c), on 
the other hand, isospin-violating New Physics effects are not parametrically 
enhanced and are generally smaller than in the Standard Model. 

For each New Physics model we have explored which region of parame- 
ter space can be probed by the — > ttK observables, and how big a de- 
parture from the Standard Model predictions one can expect under realistic 
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circumstances, taking into account all constraints on the model parameters 
implied by other processes. Table |l| summarizes our estimates of the maximal 
isospin-violating contributions to the decay amplitudes, as parameterized by 
|aNp|. They are the potentially most important source of New Physics effects 
in — > nK decays. For comparison, we recall that in the Standard Model 
a ~ 0.64. Also shown are the corresponding maximal values of the difference 
IjttK — 7|- As noted above, in models with tree-level FCNC couplings New 
Physics effects can be dramatic, whereas in supersymmetric models with R- 
parity conservation isospin-violating loop effects can be competitive with the 
Standard Model. In the case of supersymmetric models with R-parity viola- 
tion the bound ( psf ) implies interesting limits on certain combinations of the 
trilinear couplings A-^j, and A-^j,, as discussed in Ref. 127. 

7 Concluding Remarks 

We have presented an introduction to recent developments in the theory and 
phenomenology of B physics, focusing on heavy-quark symmetry, exclusive 
and inclusive weak decays of B mesons, and rare B decays that are sensitive 
to CP-violating weak phases of the Standard Model. The theoretical tools that 
allow us to perform quantitative calculations are various forms of heavy-quark 
expansions, i.e. expansions in logarithms and inverse powers of the large scale 
provided by the heavy-quark mass, nib ^ Aqcd- 

Heavy-flavor physics is a rich and diverse area of current research, which 
is characterized by a fruitful interplay between theory and experiments. This 
has led to many significant discoveries and developments. B physics has the 
potential to determine many important parameters of the electroweak theory 
and to test the Standard Model at low energies. At the same time, through 
the study of CP violation it provides a window to physics beyond the Standard 
Model. Indeed, there is a fair chance that such New Physics will first be seen 
at the B factories, before it can be explored in future collider experiments at 
the Tevatron and the Large Hadron Collider. 
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